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INTRODUCTION 


The phenomenon of reverse annealing has been observed in bulk silicon 
(ref. 1) and in silicon solar cells (refs. 2 and 3). To date, however, 
attempts to understand this detrimental effect, in terms of the defects 
formed by particulate radiation, have been unsuccessful. Recently, however, 
sensitive instrumental techniques have been used to gain increased insight 
into the behavior of defects formed in p-type silicon by particulate radia- 
tion (refs. 4-7). Hence, in the present work, we report results using rela- 
tively new information regarding defect formation in boron doped silicon 
(refs. 6 and 7) in calculations which tend to clarify and explain the re- 
verse annealing effects we observe in silicon solar cells. 


EXPERIMENTAL 


Irradiations and isochronal annealing were performed on n“^p silicon solar 
cells with 0.1 and 2 ohm-cm boron doped base resistivities. Starting mater- 
ial for the 2 ohm-cm cell was Czochralski grown while the starting material 
for the .1 fi-cm cell was vacuum float zone single crystal silicon. The cells 
were irradiated by 1 Mev electrons to a fluence of 10^^ cm""^. After irrad- 
iation, the cells were isochronally annealed in 50*^ C steps, the cells being 
held at each fixed temperature point for 20 minutes. Minority carrier dif- 
fusion lengths and AMO I-V measurements were obtained at room temperature 
before and after irradiation and after each step in the isochronal anneal. 
Diffusion lengths were measured by an X-ray excitation technique (ref. 8), 
while the AMO I-V measurements were obtained using a xenon-arc solar simulator. 
Isochronal annealing data are shown in figure 1 where it is seen that reverse 
annealing occurs for the 2 ohm-cm cell at temperatures above 150° C, but is 
absent in the 0.1 ohm-cm cell’s data. 

In order to understand these results in terms of properties of the radia- 
tion induced defects, we used a combination of our diffusion length measure- 
ments and defect data obtained from Deep Level Transient Spectroscopy (DLTS) , 
(ref. 7). The pertinent diffusion length data for the cells of figure 1 are 
shown in figure 2, where L”^ is reciprocal diffusion length squared, measured 

at room temperature after isochronal anneal at temperature T, while L”^ is 

XKR. 
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reciprocal diffusion length squared measured at room temperature after irrad- 
iation and before annealing. The DLTS, isochronal annealing data for a 2 ohm- 
cm n*^p silicon solar cell are shown in figure 3, while figure 4 shows the DLTS 
data for a 0.3 ohm-cm cell (ref. 7). In figures 3 and 4, is energy at the 
top of the valence band, while is energy at the bottom of the conduction 
band. Regarding the cells used to obtain the DLTS data (ref. 7), starting 
material for the 2 fl-cm cell was Czochalski grown boron-doped silicon, while 
that for the 0.3 Q-cm cell was vacuum float zone boron-doped silicon. 


THEORY 


Derivation of Working Equation 

To obtain an equation relating the present diffusion length data to de- 
fect parameters, such as carrier capture cross sections and relative defect 
concentrations, we use the Shockley-Read-Hall (SRH) theory of recombination 
(refs. 9 and 10). Our starting point is the SRH expression for the car- 
rier lifetime due to the i^^ single level defect. For low injection condi- 
tions we have (refs. 9 and 10) 
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Where Nj|^ is the concentration of the ith defect, n^ and pQ are equilib- 
rium electron and hole concentrations, Vjj and v„ are thermal velocities for 
electrons and holes, and are hole and electron capture cross sections 
for the ith defect and nj|_ and p^ are calculated from the relations 
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Where and Nv are the density of states for the conduction and valence 
bands k is the Boltzmann constant and Ei is the energy level of the defect. 
The thermal velocities are obtained from the relations 



Where mj^ and are hole and electron effective masses respectively. 
For P type silicon p^ is measured and n^ calculated from the relation 

tloPo = 

Where is the intrinsic carrier concentration. 

Using equations 2, 3 and 4, we find for the defects in figures 3 and 5 


^ I 


no+-Po 
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% Si:, a X 10 cm/sec ■ 

Hence equation 1 simplifies to 

Where is the relative concentration of the ith defect with 
P : .31 
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N n,a-ir being the maximum defect concentration measured by DLTS during the 
isochronal anneal. For example in figure 3, obtained from the peak 

amplitude of the + .38 defect (ref. 7). 


Since 



Capture Cross Section 

The capture cross sections obtained from DLTS (ref. 7), are listed in 
Table 1. Since the electron capture cross section for only one defect 
(E(,-0.27 eV) is available from DLTS, we cannot calculate diffusion length 
ratios from equation 8 using only this DLTS result. We have therefore ob- 
tained the necessary additional 0n values by other methods. 

With respect to the defects at Ey+ 0.38 and Ey+0.23, temperature depend- 
ent lifetime studies in 1.5 ohm-cm p-type silicon yield values for the ratio 
CJn/Op from which 0n is calculated (ref. 11). From reference 11, for the 
Ey-K).38 defect 150 and for the defect at Ey+0.23, 0n/0p ^ 1300 from 

which the values for these defects shown in Table 2 are calculated. For 
the remaining defects, values for 0^ were computed by fitting equation 8 to 
the diffusion length data of figure 2 at 250 and 300° C using Rf values from 
figure 3. In this way, the electron capture cross sections for the defects 
at E^+0.30 and E^+0.26 were obtained. Thus, the required cross sections for 
all but the defect at Ey-K).2 eV are obtained for use in equation 8. Since 
we cannot evaluate 0^^ for this latter defect, it is arbitrarily assumed that 
it can be ignored in subsequent calculations . Similar considerations apply 
to the E^+0.48 defect in 0.3 ohm-cm p-silicon (fig. 4) i.e., since we cannot 
evaluate ^n for this defect it is omitted from the calculations. 
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RESULTS 


The ratios, for the 2 Q~cm cell, calculated from equation 8, using 
relative concentrations from figure 3 and capture cross sections from Table 2 
are shown in figure 5. From the figure it is seen that the calculated and 
measured values are in reasonable agreement. From these results, the defect 
responsible for the reverse annealing in 2 ohm-cm p-type silicon is identi- 
fied as the defect at E^+0.30 eV. Additional calculations were performed 
after conceptually removing this defect, using the remaining cross sections 
of Table 1 and the remaining relative concentrations of figure 3. The result 
is shown in figure 6. 

Calculated values for the 0.3 ohm-cm cell using the data of Table 2 and 
figure 4 are shown in figure 7. These results indicate the absence of reverse 
annealing in this low resistivity cell and are in qualitative agreement with 
our data for the 0.1 ohm-cm cell. Removal of the E^+0.30 defect results in 
calculated values shown by the dotted line of figure 7. Since no DLTS data 
are available for 0.1 ohm-cm p-type silicon, we interpret these 0.3 ohm-cm 
results as an indication of trends when going to lower resistivity p-type 
silicon. 


DISCUSSION 


The good agreement between our experimental and calculated values indi- 
cates the correctness of the mathematical relationship (equation 8) derived 
from the SRH theory. Our results indicate that the defect at E^-K).30 eV is 
responsible for the observed reverse annealing. Furthermore, its removal 
could result in significant annealing at 200® C for the cells investigated. 
Since the E^+0.30 defect always appears when the E^,-0.27 defect disappears 
(ref. 7) the two defects are related in the sense that the E^,-0,27 defect 
can be called the "parent” to the 11^+0.30 defect. Hence, the atomic consti- 
tution of both defects is of importance. 

The defect at E +0.30 has been tentatively identified as a boron-oxygen- 
vacancy complex (ref. 7). The E^,-0.27 defect increases in concentration with 
increasing boron concentration and has been identified as a boron interstit- 
ial-oxygen interstitial complex (ref. 7). An alternate identification asserts 
that the E^-0.27 defect is composed of a boron interstitial-boron substitution- 
al pair (ref . 12) . A firmer identification of the E^+0.30 defect is required 

as well as additional research directed at choosing between competing identi- 
fications for the E^-0.27 parent defect. Such identifications would be an 
invaluable guide to processing efforts aimed at decreasing the concentration 
of these radiation induced defects. 
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TABLE I.- CAPTURE CROSS SECTIONS FROM DLTS 


DEFECT ENERGY 
LEVEL 
(eV) 

CAPTURE CROSS SECTIONS 

(cm2) 

ap 

ON 


2x10-16 

- 

Ec-0.27 

- 

2x10-13 

Ey+0.23 


— 

Ev+0.30 

2x10-16 

- 

Ey+0,26 

__ 


Ey+0.2 

- 

— 


TABLE II. - ELECTRON CAPTURE CROSS SECTIONS 


ENERGY LEVEL 
OF DEFECT 
eV 

ELECTRON 
CAPTURE CROSS 
section (cm2) 

METHOD FOR 
OBTAINING CROSS 
SECTION 

E^-0.27 

2x10-13 

DLTS 

Ey+0,38 

3xlO-W 

•an/op 'V* 150 

E/0.23 

3.9x10-13 

*aN/op 'v* 1300 

Ev+0.30 

3.6x10-13 

FIT OF Jy RATIO TO DLTS (25CPC) 

E^+0.26 

g.Sxio-w 

FIT OF ^ RATIO TO DLTS (30CPC) 


•COMPUTED FROM cN/ap RATIOS FOR 1.5 n-cM P-TYPE (FZ) FROM SROUR, 
CURTIS. OTHMER. CHIU & DEOKAR: REPORT HDL-TR-171-^-U.S, ARMY 
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Figure 1.- Isochronal anneal of n'''p silicon solar cells after 1-MeV-eiectron irradiation, 
time at temperature. 20 min. 
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Figure 2.- 1/L‘ ratios from diffusion lengths measured during isochronai anneal. 
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Figure 4.- DLTS data for 0. 3ohm-cm P-SL 
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Figure 5.- Calculated and measured 1/L^ ratios for 2 ohm-cm cells. 
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Figure 6.“ Calculated and measured 1/L^ ratios, omitting Ey + 0.30 effect. 
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O 0.3 -CM CALCULATED 



Figure 7. - Calculated and measured lIL^ ratios for low-resistivity P-Si. 



